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IMPACT THEORY FOE SEAPLANE LANDINGS 
By Stanley U» Benscoter 

SUMMARY 


A step landing of a seaplane on smooth water is analyzed by 
solving Newton's differential equation of motion. Various limiting 
assumptions are made for convenience in the analysis vtfien they 
appear reasonable. Formulas are developed for the maximum accelera- 
tion, the maximum draft, and the draft, at the' instant of maximum 
acceleration. The solution is dependent upon the definition of 
the apparent -water mass and indicates the need for further basic 
research in regard to this water mass. Graphical representations 
of approximate solutions for design, use are presented. Comparisons 
of theory with experimental values obtained in the National Advisory 
Committee for Aeronautics impact basin are also illustrated. . 


U3TB0DUCTI0N 


In the analysis to be presented it will oe assumed that the 
seaplane is landing on smooth water. If the reader has ever 
"skipped" a flat rock off the surface of a smooth lake or pond, 
he has witnessed a physical action which is similar to that of a 
seaplane landing. If the rock is thrown horizontally, close to 
the water surface, it will "skip, “ or rebound, frcm the water 
surface several times before final submerging. Similarly, in a 
seaplane landing, if the trim were held constant, several rebounds, 
or impact periods would occur before the seaplane finally settled 
onto the water. During each of these impact periods the reaction 
force exerted by the water on the hull increases from zero to a 
maximum and then decreases again to zero. In each consecutive 
impact period the maximum value of the reaction which is developed 
is less than in the previous period. Thus the maximum value of' 
the reaction force will occur during the first Impact period. 
Consequently a mathematical analysis may be limited to a study 
of the motion and associated forces during .the first impact. 

During this period the lift force on the wing Is approximately 
equal to the weight of the seaplane . As a matter, of. . convenience, 
it will be assumed throughout the analysis that the wing lift and 
seaplane weight remain equal through -the impact period. 
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• SYMBOLS- 
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m 

r 

0 
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V 

x 

z 

a 

3 

7 

€ 

X 

3 

P 
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loaded area 
lift coefficient 

half vidth of a unit loaded atrip 
total reaction force 
force oh a unit strip of loaded- area . 
acceleration due to gravity - 

length of loaded area 
mass of seaplane 
additional mass of water 
* ratio of tan 7 tc tan t 
step draft normal to keel 
time 

velocity - - 

displacement parallel to keel 

displacement normal to keel 

factor occurring in the definition of € 

dead-rise angle 

flight- path angle 

coefficient which defines a 

length parameter ' equal to '\~f M/p 
‘ ratio of m to' M - 
mass density of ' -water - 

' *» .'*■ 

trim a^agle ' 

convenient universal function 
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c- .space coordinate parallel to the X-axis 

t space coordinate parallel to the Z-axis 

Subscripts: - 

h;.*“ •' ' horizontal component ' 

f . t 

v Vertical component 

' m value of quantity at instant of maximum acceleration 

n value of quantity at instant of maximum draft 

o value of quantity at instant of entry (t, a 0) 

Prime i A prime is used to indicate differentiation with respect 
to a space variable. 

Dot: A dot over a letter is used to indicate differentiation 

with respect to time. 

Bar: A bar over a letter is used to indicate that the quantity 

is of two-dimensional nature. 

DISCUSSION OF MOTION ' : , • 


At the instant of entry of the seaplane into the water there 
are certain initial conditions of the motion. These are illustrated 
in figure 1(a) . The trim angle T, made by the keel with the water 
surface, is assumed to. remain constant during the impact;. " "This" 
assumption is approximately true in a step landing. The initial 
velocity is V Q at the instant of entry and the direction of ■ V Q 
is given by the initial flight -path angle y 0 . The enlarged view 
of figure 1(b) shows the seaplane entering the water at point 0 
on the water surface. The point 0 is to be chosen as the origin 
of a fixed reference coordinate system. The point A on the 
seaplane is defined as the intersection of the keel and the step. 

At the instant of entry (t = 0) the point A of the seaplane 
enters the water at point 0. The motion of .point A may be 
regarded as the motion of the seaplane.. 

In figure 2. (a) the shape of the path of motion.- -of point A 
is indicated. The motion of point A along this path will be 
referred to as the "motion of the seaplane throughout' the analysis. 
The seaplane enters the water at point 0 Snd emerges at point P. 
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At some intermediate position the draft passes through a maximum. 

If the differential equation of motion of the seaplane can be 
properly developed and Solved, the relationships between time, 
displacements, velocities, and accelerations during the Impact 
can be determined. The bottom reaction will also vary during the 
motion somewhat as illustrated in figure 2(b) . This reaction is 
proportional to the vertical acceleration of the seaplane when its 
wing lift Is equal to its weight. It Is found, both from theory 
and experiment, that the maximum value of the upward acceleration 
(maximum deceleration) occurs before the Instant of maximum draft. 
In experimental research work oh float models it is possible to 
measure the maximum draft and the draft at the instant of maximum 
acceleration,- as well' as' the maximum acceleration. A correct 
theory should predict all three of these quantities in agreement 
with experiment . - 

- The motion will not only be affected by the initial condi- 
tions but also by certain physical properties of the seaplane. 

The properties which are of primary Importance are the mass of the 
seaplane and the ’shape of the hull bottom. The elasticity of the 
seaplane structure is of secondary importance in most cases and 
will be neglected in this analysis. In the special case of a 
large flying boat',’ with engines in the wing, the elasticity of 
the wing may have an appreciable effect upon the motion. The 
bottom of the hull is of wedge shape with some curvature, or 
flare, of the bottom. As a matter of convenience it is assumed 
In the analysis that the bottom is a flat-sided wedge without 
flare. The bottom shape Is then defined by the dead-rise angle. 

The two physical properties of the seaplane which affect the motion 
are, then, the mass M and the dead-rise angle. 3- 

' In a landing with or without power the Initial conditions of 
the motion may vary over a range of values depending upon the 
pilot's skill and landing technique. Consideration of reasonable 
landings, for design purpose's would require a considerable amount of 
statistical evidence regarding -the Initial conditions of motion 
as obtained In actual flight operations. Such evidence is not 
available at present. Plight research in this field is seriously 
needed. > 

If -the friction of the water on the hull bottom is neglected, 
the component of. velocity of the seaplane parallel to the keel 
will remain constant during the impact period. In experimental 
research work on model floats the float is launched from a heavy 
carriage which rolls 'on a horizontal track, usually overhead. 

This method of launching the' model allows it to have freedom of 
vertical motion but forces It to maintain almost a constant 
horizontal velocity because of the- large inertia of the carriage to 
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which it remains attached. In the present analysis it will he 
assumed that the horizontal component of velocity remains constant 
during the impact in order that the. resulting formulas can he 
compared directly with experimental evidence. Since the trim angle 
in a step - 'landing must he small, the maximum reaction given hy a 
solution involving constant horizontal velocity -will differ only 
slightly from a solution for constant velocity along the keel and, 
hence, should he satisfactory for design purposes . 


DISCUSSION OF REACTION DISTRIBUTION 


When a hull is immersed a small amount at a particular instant 
during the impact, the loaded area is of approximately triangular 
shape. The "base of the triangle is at the step. Due to the effect 
of "piled^up" water the triangular loaded area is larger than 
would he determined from the geometrical intersection of the hull 
with the plane of the original water surface. For design purposes 
it is necessary to know the transverse distribution of pressure on 
this area and the longitudinal, or fore-and-aft, distribution. 

The longitudinal distribution must he rounded off in some manner 
at the step as indicated in figure 3(a) . This three-dimensional 
effect can only he brought into the analysis, at present, hy means 
of an aspect-ratio factor applied to the pressure at all points 
of the loaded area. 

In figure 3(h) a unit transverse strip of the loaded area is 
shown with the load acting on this strip. The reaction force on a 
unit strip at a particular station is dependent upon the draft 
of the keel at this station as well as the components of velocity 
and acceleration of the seaplane normal to the keel. The component 
of velocity parallel to the keel does not affect the force on :&he 
unit strip. In order to determine the total reaction, the reaction 
per unit of length must he determined and then an integration 
performed in the fore-and-aft direction over the loaded area. The 
reaction force on a unit strip can he determined hy considering 
the two -dimensional case. For this reason the two-dimensional 
case will he analyzed first. This analysis will he carried out 
more completely than is necessary for the purposes of this paper 
in order to provide formulas which may he useful for comparison 
with results from future experimental research studies.. The 
transverse pressure distribution on a unit strip will he somewhat 
as shown in figure 3(c) • No attempt will he made in this paper 
to determine' this distribution. 
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TWO -DIMENSIONAL CASE 

- ‘ ■ *\ * 


In this case a wedge -shaped' float’ is assumed to enter the 
i3ater' : vertically as shown in figure k. The wedge is assumed to 
he of infinite length and infinite width. The action of a unit 
Slice of the wedge and a unit slice of the fluid is considered iji 
the analysis. -It- is assumed that the momentum of the fluid is 
equal to the product of an "additional mass" St and the 
instantaneous velocity z of the float. The displacements, 
velocities, and- accelerations are considered positive downward.' 
The additional mass is a half cylinder of water as shown in 
figure h. This elementary concept of the additional mass and its 
use ‘in the calculation of the reaction on the float were firBt 
given hy Th. von Karraan (reference 1) . The assumed additional 
mass is in agreement with the additional mass which is known 
(reference 2) to he associated with a flat plate moving through 
an ideal fluid of infinite extent normal to the plane of the . 
plate. . 

Due to the effect of piled-up water which would actually 
occur and the fact that the float has a wedge "bottom rather than 
a flat "bottom, it is necessary to multiply the additional water 
mass as -shown in figure h by correction factors.- Such factors' 
were first introduced by H. Wagner (reference 3) . In any case, 
however, the water mass Hr must be proportional to the square of 
the displacement and may be written as 

5 =* ~i z 2 : (1) 


The bar over a quantity is used to Indicate that_the quantity is 
of two-dimensional nature. The derivatives of m with respect 
to displacement and time are required in the analysis and may be 
written in- the following forms : 


dm 

dz 


m’ 




( 2 ) 


dm 

dt 


5 


=> 2ezz = 


2m& — . 


= 


(3) 


The reaction force is due • to the rate of change of momentum 
of the fluid. The change of momentum is, due both to the change of 
velocity of the float and the change of size of the additional 
mass caused by a change of the loaded width. All effects of 
viscosity and buoyancy are omitted from the analysis. The force 
may be expressed in the following form: 
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f = ~ (mz) = mz + mz ( 4 ) 

' l 1 •’ r ~ "■ . ’..T i * 

Using equation ( 3 ) this may "be written/' ” 

• v • .. <• f >■ mz,+ m'fc 2 ’ • ; (5) 

This force acts upward and, thus it must he preceded by a 
negative sign in- ■writing Newton’s law of action for -the. float, 
Newton’ s law, when applied to the. float, is expressed hy, ’ ■ 

m m -f ( 6 ) 


Tfya solution of this equation will he that which was originally 
given hy von Earmsfcn. Substitution of equation (5} gives, 



+ mir 


* p 

+ m «* 


( 7 ) 


(M + m) z 


.2 

-+ m z 




It becomes convenient in the ana.lysis to introduce the 
quantity n as the ratio of the water mass to the float mass. 


V- 



( 9 ) 


Equation ( 8 ) becomes. 


(l + JT)z + n'z 2 = 0 


( 10 ) 


From equations (3) 

and (9) it Is seen that, : -- 

** 


{!’ z . 

(11) 

Hence equation (ID) 

may he written, 



dz ■ . djl 

(1 + f) 2t + z dt“° 

“ (12) 

'or 

... •. 



(1 + fT) dz + £ djl = 0 

( 13 ) 


"■ ■ ' ■' ‘ ■ - : ■ 

■ 

* 

dz . dil 

z 1 + il ■ "" 

• , ' . m 
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Integration gives,. 


log z + log (1 .*»?). =» ; log;.C . ■ .QL5) 


The constant of integration is determined from the Initial con- 
ditions. ‘When t s* 0 the displacement z is zero. Hence the 
mass m is zero and also the ratio |i. The initial velocity at 
the instant of entry may be indicated as ‘£ 0 . Then ■ C «.£ 0 ' • 

and equation '41-5} ‘"becomes. ' " ' * 

, 

log z + log (1 + jT) =» log & Q (16) 


or 


o 

(T+7tT 



This equation vas first given by von Kerman (reference 1) . TMb 
differential equation may be regarded as a formula for the velocity 
in terms of the displacement. Correspondingly, equation (1C)‘ may 
be regarded as a formula for the acceleration in terms of the 
displacement and velocity. Solving equation (10) for the accelera- 
tion gives, .... 


-U*z 2 ’’ 

(1 + ft 


(18) 


From equations (2) and (9) 


it i$ seen that. 



Equation (l8) becomes. 


(19) 


z 


-2jlz 2 ■ 
z(l + p) 


( 20 ) 


Equation (17) may be substituted to obtain the acceleration as a 
function of the displacement 


-2UZn 2 

z(l +_ il)3 


( 21 ) 


This formula -was also given by von Harman (reference 1) . By 
assuming various values of z the velocities and accelerations 
can be computed from equations (lj) and (21) . 
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In order to illustrate the variation of velocity and accelera- 
tion -with displacement a simple numerical example has "been calculated 
for the following physical constants. 

7 * 0.05 slug per cubic inch 

M =» 1 slug per inch 

z 0 .,“ 6 feet per second . 

■ ■ . . - . . , _ . i 

A plot of i and & against z is shown in figure % It is 
seen that the acceleration passes through a maximum at a displace- 
ment of 2 inches. The accelerations are plotted in terms of the 
acceleration g due to gravity and. the ordinates to this curve 
may "be multiplied "by the weight" of the float to obtain the reaction 
force. It may be noted that the velocity remains positive 
indefinitely. This means that the float moves continuously down- 
ward without reversing its direction of motion. This is due to the 
omission of the buoyancy force from the equation of motion. At the 
displacement for maximum acceleration the buoyancy force is very 
small and, hence, it appears reasonable to emit the buoyancy force 
in determining maximum acceleration. 


For purposes of "structural design it is not necessary to . 
know the value of the acceleration at various instants but only 
the maximum value. This may be obtained by differentiating the 
formula for the acceleration. Differentiating equation (21) with 
respect to z gives. 


. t (x . 

4a a \ 1+iy 


( 22 ) 


Setting the quantity within the brace equal to zero and solving 
for (A gives. 



(S3) 


The subscript m is used to indicate the value of a quantity 
when the acceleration is a maximum. A universal constant value 
of 1/5 is obtained for {r m which holds true for all values of 
M and (3. This constant was first published by J. Sydcrw 
(reference 4) . 

From equations (1) and ( 9 ) it is seen that the displacement 
at maximum acceleration is. given by. 
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or 



Substitution of equation (23) into equation (21) gives the formula 
for maximum acceleration. 



Equations (25) and (26) give tjie. maximum acceleration, in terms of 
the physical properties of the float and, the initial sinking speed. 

In order to plan proper instrumentation in experimental 
research work it is desirable to know the time which elapses from 
the instant of entry to the instant of maximum acceleration. The 
relationship of displacement, velocity, and acceleration to time 
may be determined by integrating equation (17) . 



(1 + jl)z « z Q 

(27) 

or 




(1 + -) iz - i a dt 

(28) 


Direct integration gives. 



"ez^ t , 

z + $r " z o t 

(29) 

or 


(30) 

The constant, of integration is zero, The time 
tion is obtained by substituting equation (23) 

for maximum accelera 


■16 / z m\ 

55 In) 

(31) 
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Equation C 30) is a cutic equation in z which could be solved 
for z in terms of t. However, -it. is more convenient for calcula- 
tion purposes to assume various values of . z and compute the 
associated values of'-, t-. The associated, values of t and z 
can also he computed. Displacement, velocity, end acceleration 
can then he plotted against time. Graphs of this nature have 
been -calculated for the numerical example previously considered 
and are shown in figure. j(?. • : - r. • , 

•The two-dimensional case has been given considerable, attention 
because the steps which have been taken suggest similar steps 
which can be taken in the three-dimensional case. 

THREE -DIMENSIONAL CASE 


The ‘three-dimensional analysis will be made for an. ideal . 
prismatic float of semi -infinite length and infinite width with 
a flat-sided wedge bottom. The float begins at the step and 
extends indefinitely forward. It is assumed to have no. afterbody. 
The float is assumed to maintain a constant trim angle throughout 
the impact period. 

In figure 7 an- ideal float is shown immersed in the water 
after having travelled along a path of motion indicated by the 
dotted lin e. The origin of coordinates is chosen as the point 
on the water surface where the point, A of the float enters the 
water. The X-axis is chosen to be parallel to the keel and the 
Z-axis is normal to the keel. The positive directions of the 
axes are as shown. In the position of the float as shown the 
flight-path angle has .decreased, from the initial value. 7 0 . to 
the Instantaneous value 7 . 

The fluid beneath the float is shown divided into unit 
slices by planes which are perpendicular, to the keel. -Due to ■■ 
the lack of a satisfactory three-dimensional hydrodynamic theory- 
it is necessary to assume that the reaction on the float is due 
to the motion in these slices of fluid which have a combined 
width equal to the wetted keel length. It is necessary, to assume 
that the motion of the fluid in each slice is independent of the 
motion in adjacent slices and that, the displacements, velocities, 
and accelerations in each slice are in a plane normal to the keel. 

A particular unit slice has been indicated by cross batching 
in figure 7 • At the center line of this slice the keel is at a 
distance . £ -from the water .-surface. ' Corresponding to the two- . 

■ dimensional theory previously ’developed there may be associated ■ 
with a particular slice of fluid the additional water mass 5 
given by 
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5 = e £ 2 . 


(32) 


It' ig understood that, -the coefficient ' 7. cow contains- a Ire duct ion 
factor for aspect-ratio effect. 

The center line of the cross -hatched- sli ce .in figure 7 is at 
a distance £> . - from the leading 'edge Of the loaded area. The 
quantities g and £ are instantaneous values. The fluid in 
the cross-hatched- slice exerts, on the float, a force f given by 
the f allowing fcnmula: • ' r - 


f n 


d_ >- 

at 


(m£) « + m '£ 2 


(33) 


In order to obtain the total force on the entire bottom 
of the float, the above formula may be integrated over the wetted 
keel length I*,.. . " ■ .... 


F 


fit r>L .. _ LO 

f d£ : ». (m£ + m*t 2 ) a£ , ,(34) 

0 . Jo- 


Substitution of equation ( 32 ) and its derivative gives, 

+ jo 


F=*J d§ ■ ^ * i2 


(35) 


The distances £ and g are related through, the trim angle. 

. £ » £, tan t ( 36 ) 

Substituting this formula into equation ( 35 ) and bringing the 
constant terms outside of the integrals giyes^ 


■ pL 

F = i"£tan 2 T I <| 2 d£-+ 2c t 2 tanT / 

Jo Jd 


et tan 2 T ' 


+ 27 £ 2 tan t | 


(37) 


(38) 


■ ,,-f 


It may be qeen from the definitions of £ and z that, although 
their magnitudes are different, their first and second derivatives 
with respect to time are equal. * 


t 


£ 


( 39 ) 
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Substitution in equation ( 38 ) gives, 

F = a * (et,2 tan t). i 2 


(40) 


If s is the value of £ at the step, it. will he - related to the 
length ■ L ' through the trim angle . 


■ s = .-Jj- tan T 

Thus equation (40) may he written, 


F 




2 + (£Ls )£ 2 


(41) 


(42) 


The quantity which is within the first parenthesis may he recognized 
as the mass of a half cone of water. -This three-dimensional addi- 
tional water mass is illustrated in figure 8. This mass may he 
written as. 


m =» 


St>s 2 


.3 tan T j 


s J •* e S' 


(43) 


The derivative of m with respect to s may he written, 

• r_2 ■ 


dm » ‘ €s“ 

ds * m tan r 


3es' 


.2 


(44) 


The q uan tity which is contained in the second parenthesis of 
equation (42) may he recognized as - equal to the above derivative. 
Hence equation (42) may he written, 


F = m2 + m*2 2 


(45) 


This fo rmula was originally, developed . in approximately the same 
form by W. L. Mayo (reference 5) »- , 

In the above equation it is of particular importance to 
note that the derivative m* is taken with respect to s rather 
than z as in the two-dimensional case. The importance of this 
dist inc tion arises from certain new features of the three- 
dimensional case. Although -the reaction of a particular unit 
slice on the float is not affected hy the component of velocity 
parallel to the keel, the number of unit slices and the rate of 
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change of tills number are affected. "by this component of velocity. 
Or, from another 'viewpoint, the additional water mass is determined 
hy the step distance s which varies due to & as well as z. 

This matter has "been given rather extensive discussion by Mayo 
(reference 5) * 

The analysis will be based on the assumption that the hori- 
zontal velocity Tij remains constant. This means that the 
resultant acceleration is vertical. Hence the equation of motion 
must be written for forces in the vertical direction. The 
resultant value of the acceleration, which is vertical, is given 
by t/coB T. The vertical component of reaction force is given 
by F cos T. This assumes that the resultant reaction force is 
normal to the keel. Newton's equation of motion may be -written. 


Mz 

COS T 


+ F cos T 


0 


(46) 


Substitution of equation (45) gives, 




COS T 


+ mz cos T 


p 

+ m'z 


COS T 


(47) 


An equation of this type, showing considerable variation, however, 
was originally published by Mayo (reference 6) . A single integra- 
tion of .his equation was performed by Mayo after which relations 
between displacement, velocity, acceleration, and time were 
obtained by numerical methods. 

It is convenient to introduce the dimensionless constant i-i as. 


2 

m cos T 

* = —w— 


(48) 


Equation (4?) may be written as 


or 


t- ( m cos 2 t\ /m*. coe 2 r\ ±2 ^ 

fe;+ \ ; — ' ir~7 z + — / 2 * 


(i + \i)z + z 2 ^ = o 


O' • (49) 


(50) 
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This equation nay he ^rritten in the form, 

• ■' . ■ * .* - r - . ■ — J- \ 

- - *“ ’ r ~ 

/' .._-J 1 + ^ )2+ .(4)«“ 0 (51) 

It is now necessary to develop expressions for i and s in 
terms of the constant velocity 7^ and a convenient variable. 
The following geometric relations are apparent from figure 8. 

z ■ 7 sin (7 + t) (52) 



« 7 sin 7 cos T + 7 cos 7 Bin T 


7 V cos T + 7^ sin T 
■(7 t qos t + 7 h sin t) 


2 


fv cost Y 
+ x ) 


ein^T 


(53) 

/ • 

( 5 *) 

(55) 

(56) 


The ratio 7^/V-^ is the tangent of the flight-path angle. It is 
convenient to introduce the dimensionless, variable, r . given by 
the ratio of tan 7 to : tan t . 


tan 7 
r ~ tan T 


(57) 


■'The flight -path angle may be regarded as a dy n a m ic trim angle and 
"'the ratio ■ r is the ratio of the dynamic to the static trim angle 
•' "Eq.uata.qir (56) .becomes^ ' 

(z) 2 => (r + l) 2 7jj® 'sin^T ' : " (58) 

The depth of the point A oh the float below the water surface 
at any instant is equal' to'. . s cos' T. The vertical velocity of A 
may be obtained by differentiating this quantity with respect to 
time . Hence, 


7 V = s cos t • ; ( 59 ) 




1 6 


KACA TN No. 2>37 


It is desirable to determine the derivative of r with' respect 
to time . 



(62) 

(63) 


Substituting equations (58), (6l), and (63). into equation (51) 
gives, 

(1 + u) V h sin T v h sln T C os s T ^ * 0 (64) 

or 

(1 + |i) dr + gos 2 t dp. ■* 0 (65) 

r 

Separating the variables and integrating gives, 

log (1 + r) + 3~7"'r + cos ^ T log (1 + u) **0 (66) 


The constant of integration C must be determined from the initial 
conditions. At t =» 0^ the distance s is zero, the water mass m 
is zero and, hence, (J. «* 0. The value of r at t = 0 may be 
indicated as r ? . This value will be determined by the initial 
flight-path angle 7 0 • 



tdn Tq ■ • 1 -v 

a ■- 

tan T 


(67) 


Substituting the initial conditions in equation (66) gives the 
value of C. 

C » log. (1 4> r D ) + (68) 

’ x + r o 

Equation (66) becomes, ‘ " . ' 

log (1 + r) + 577-7 + cosS- log (1 + fi) *». log (1 + r 0 ) + 
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This 'equation expresses a relation between instantaneous values 
of u and r for a particular trim angle and a given initial 
flight-path angle. By assuming various values of r from r 0 
to zero the associated values of p, can readily he computed. 

This provides a direct relationship "between the vertical velocity 
and the draft of the float. 


At this point in the analysis it is necessary to return to 
the original differential equation of motion, just as in the 
two-dimensional case, and to regard- this equation as a formula 
for the acceleration in terms of the displacement and velocity. 
Solving equation (50) for 2 gives, 


z g dp 
1 + p ds 


(70) 


From equations (44) and (48) 


it is seen' that. 


ds • s 


(71) 


Substitution of equations .( 58 ) and (71) into equation (70) gives, 

£ - 1 ( l 2 ^) (1 + r)2 \ 2 3tr?T 


A relationship between p . and r ■ which holds true at the instant 
of maccimum acceleration can be obtained by differentiating the 
formula for z with respect to s and equating the result to. 
zero. Ih order to carry out this differentiation it is necessary 
to establish the derivative of r with respect to s. Dividing 
equation ( 63 ) by : -s gives, : 


dr dt z 

dt ds s Tjj sinr 


(73) 


Substituting equation ( 6 1) on the right-hand side gives 


dr a z / cos g T \ = H 
ds T h - sin T i B ^ n 7 , rV^ tan g T 


(74) 
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Differentiating equation (72) , and using the above formula, results 
in the following equation, 


48 z 

2 - 4 ' 

. 64 . (l + r\ 0 

33 "** 

da s 

1 + 4 

1 + 4 \ r } 


(75) 


The quantity within' the brace may bo equated to zero to obtain a 
relation between 4 and r which must hold true at the instant 
of maximum acceleration. These particular values may be indicated 
as Urn and r m . 


? - 4. 


m 


64 


m 


r l + r n 


1 + bn 1 + bn V r ) 


cos^r 


m 


(76) 


This equation may be solved for 4 m or r m to obtain the formulas 




2r. 


bn “ 


m 


r m (l + 6 cos^r) + 6 cob^t 


(77) 


64 


r =» 
m ■ 


m 


o 

COS T 


2 - (1+6 cos {? t ) 4, 


m 


(78) 


By integrating the original differential equation of motion 
a relation between 4 and r was established as given by equa- 
tion (69) • This relation must hold true at all Instants of the 
impact period. Equation ( 76 ) gives a relation which must hold 
true at the Instant of maximum acceleration. Equations ( 69 ) 

■ and ( 76 ) may be solved simultaneously to obtain solutions for 4 m 
and r m in terms of the known initial quantity r 0 . 'The resulting 
values may then be substituted into equation (72) to determine 
the maximum acceleration* 


Since equation ( 69 ) is of a transcendental nature, no direct 
analytical solution can be obtained for the maximum acceleration. 
If B m and Zjjj are introduced to Indicate the values of these 
quantities when the acceleration is a maximum, equation ( 72 ) gives. 


m 


- 1 / \ \ 

% [} .+ bn / 


(1 + r m ) 2 sinST 


(79) 



NAGA TN Ho. 1437 


From equations (43) and (48) the distance is given by. 


m . \/ £ COS c T 

f 

It is now apparent that, "by using numerical and graphics.! methods 
for determining ■ u m and r a , it would he possible to calculate 
the maximum acceleration if the parameter e were Isnown. This 
parameter determines the apparent mass of the water. No attempt 
will he made in this paper to discuss the various theoretical and 
experimental studies which have been made to determine this 
parameter. These developments have been discussed by Mayo (refer- 
ence 5) • A considerable amount of basic research work in this 
field is still needed. 

The maximum acceleration may be expressed in terms of various 
components of the initial velocity. That component which seems 
most rational to use is the component normal to the keel; which 
is indicated as i Q * From equation (5®) it is seen that, 

£ 2 

V - g ( 81 ) 

(1 + r D )^ 

Substitution into equation (79) gives. 


t(3 


3- + 

1 + r. 


For design purposes it is most convenient to express the 
acceleration in terms of the resultant initial velocity V 0 . 
Using the definition of r as given by equation (57) it is 
readily shown that, 


p p 

1 + r Q taut 


Substitution into equation (79) gives, 

g m .1 / ^m V (U r n> 

Zm Sm I 1 + »W 'C 1 + r c 


2 V 0 2 Blxfr 


ta n 2 t) 
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■ At' ' th& instant of - milinum 4 xe.it ■ 'the flight -path • angle - is" zero . 
Substituting r ^0 in equation (69) gives, 

’ f ■. * •£ 

cos 2 t log (1 + 1%) '*> log (1 + r 0 ) - - — 2 — (85) 

_ 1 + r o 

' ► - ... ,* , „ . t . _ , . - - 

• The subscript -n is used to indicate the value of functions at 
'the instant of maximum draft.. For any given value of • r Q , the 
corresponding value' of p n can be computed. The- maximum draft 

• may .then be determined from equation (80) by changing the subscript 
m to- 


approximate" formulas for design 


■ Since the trim angle is- a smell angle in normal landings it 
appears reayo-urble to assume: cos- T » 1 wherever this proves to 
be convenient. Under this assumption- the approximate form of 
equations ( 6*9) and (77) becomes. 


log (1 + r) + + log (1 + M-) » log (1 + r 0 ) + 


1 + r_ 


( 86 ) 


2r 

m 


m 


7 r m- + 


(87) 


A special case is worthy of note for research purposes. In the 
case of a vertical drop, with guides to maintain ■=> 0, the 
value of r will be infinite at all Instants, Dividing 
numerator and denominator of equation (87) by r m and setting 
IV ® a. gives the universal constant u m » 2/7 for the three- 7 
dimensional case. This constant was first arrived at by 
I. W. McCaig of Great Britain. 


Equations (86) afid (87) must be solved simultaneously for 
various values of r Q i A graphical method has been used for this 
purpose. An illustration of a particular solution is given in 
figure 9 . Curve number! is a graphical representation of 
equation (86) for r_ * 1 . Curve number 2 is a representation 
of equation (87) • ■ The intersection of the two curves gives the 
' values of and r_ -for r p = 1. I-n order -to show the small 

error involved in neglecting the' effect of trim angle, the 
dotted curves have been plotted in figure 9 from equations (69) 
and (77) for t «* 8°. 
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Using the method illustrated in figure 9 , the functions u m 
end r a were computed to give the graphs shown in figure 10 . 

These universal curves are independent of the physical properties 
of the float or the coefficient which defines the additional water 
mass, 

The ^ply faotor remaining to he determined in the formula 
for acceleration is the step draft Sq. This quantity is affected 
hy the float properties and the water mass coefficient. The 
additional water mass may he written in the following form: 

a * £s 3 m (a^ 038)^0 (88) 

The factor is introduced to account for the primary effects 

of t and 0 upon m. This factor is determined from tfee inter- 
section of the float "bottom with the original water surface. It 
defines a half cone of water which is the three-dimensional analog 
of the original two-dimensional concept of von Kdrm&n. Thus, 
from the geometry of figure 8 it is seen that, 



7f- 

6 tan t tan 2 0 


( 89 ) 


The factor otg^ i fl introduced into equation (88) to account 
for secondary effects of T and 0 upon m. Some of these effects 
ares 


(1) Effect of piled-up water. 

(2) Dead-rise angle correction to flat-plate theory, 

(3) Aspect -ratio effect. 

There is a serious need, for both theoretical and experimental 
research s'tudies to establish these factors accurately. Using 
the best evidence available, Mayo (reference 5 ) has suggested the 
following formula j 

^3 . o.te Wf> - if ( 1 -jSSUy) .(90) 

Eor a Practical range of t from 3 0 to 12° and 0 from 20° 
to 25 , equation (90) gives 0^3 the narnga . 0.9 < a ^3 < 1,3. 
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' la: order to obtain dimens ionlee a quantities for plotting 
purp&Ses, it is convenient to introduce the fundamental length 
parameter . ■ 



(91) 


Dividing through equation- (88) by M, and usings equation (91) 
gives the dimensionless’ variable p. as, 

" ■ ' 11 *( a l ot B ?) 3 ( >" ) 


or 




By using figure 10 for p the right-hand side of this equation 

: has been computed for 0,3- 0.9, 1.1, and 1.3. These theoretical 
curves are shown in figure 11. For a given float, landing at a 
given trim angle and flight-path angle, the step draft s m may be 
•.determined to a practical degree of accuracy from the solid line of 
figure 11; " • • .. ...... 

For plotting purposes, equation \ 84) ’’ may - be written in the 
following form by using equation (93) • * 


M, 

^7 


<*2 / 3 % \ (1 + r m )^ b ± 7 t 

yffia \ k 1 + ti m) (1-f r^.tan^-r) 


(94) 


The- right-hand side of this equation is a function of r D which 
may. be- computed and plotted for various, trim angles. Using an 
^average value of a 1.1 this function is plotted in figure 12 
for t = 3° and in figure 13 for t = 12°.' Frcp such curves the 
maximum acceleration may be readily determined for- any float. 

The reaction force is then obtained by multiplying - the acceleration 
by the mass of the float. 


The longitudinal distribution of the reaction may be determined. 
The .ordinate to .-the load distribution diagram is given by f. 
Referring tc^'equation (33) it is seep that ,_f > consists of two p 
parts. The' first "part is proportional’ to m which ..varies as & . 
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Hence the first part -has a quadratic variation. The second part 
has a linear Variation since .S’ varies as £. ■ Integration of 
the first part of f - gives the- first part or the total reaction S’. 
The magnitude ,of this first part of "7 - Is' seen to te " mfi or 
Consequently the second part of 3? taist he given by - -(l + u)M£; 
The two parts of the reaction and their distribution are indicated 
by the diagrams" of figure It. •' = -- 


COMPARISON 07 THEORI WITH EXITSIMENT 


Some excellent experimental results have been recently- 
published by the National Advisory Committee for Aeronautics for 
landing impacts on smooth water These teste were run with a 
float weighing 1100 pounds and having a dead-rise angle of 22 . 5 °. 

The floats were launched at various trim angles and various 
flight-path angles. Experimental values will be taken from two 
reports by S. A. Batters c® (references 7 and 8) for T *» 3° . 

and. T. * 12°. No attempt will be made to explain the discrepancies’ 
between theory and experiment, -These discrepancies may be due to 
any of .a number of causes which ’are not sufficiently well under- 
stood at- present. * 

In figure 11 experimental values of the relative step draft 
at maximum acceleration have been plotted. The scattering of the 
points indicates the necessity for accurate ■ displacement measure- 
ments and the need for effective statistical methods of interpreting 
hydrodynamic test data. The experimental spread of the points is 
greater than the theoretical spread predicted by the Mayo formula 
for og. The general trend of the points for T = 3° is higher' 
than that for t » 12°. This is in direct opposition to theoretical 
prediction and there is no apparent explanation of the contradiction. 

In figures 12 and 13 experimental values of relative accelera- 
tion have been plotted. The agreement between theory and experiment 
is very S°od. The curves are plotted for an average value of ctg. 

The value of at the instant of, maximum 'draft may be 
computed from the following approximate form of equation ( 85 ) . 

- v ' r 

- ■ . . log (1 + hn) ' * log .(1 + r o ) - ( 95 ) ’ 

The right-hand- side of- equation {93) has been computed for three 
values of ctg to give the curves shown in figure 15 . Experimental 
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values of the . relative draft ;are .also 'shown. ’ In general, the 
theory predicts.larger draft values, .than are found from- expert - 
ment . The same contradiction "between theory and experiment is 
found for the maximao draft values '-of figure 15 as was J noted in 
connection with .the draft values shown in figure Hr- ' ' 

I".-" "• • "• ' ' ' " ' . ‘" 

Corre spending to equation (94) other- dimensionlesa accelera- 
tion formulas may "be derived from equations (79) and (82) , 



It. may "be noticed that the right-hand side, of each of these equations 
is affected by variations .in t . and . p only through the' minor 
variations of ' the secondary factor ctg. Using an average value 
of ot2, a theoretical curve has been computed for each of "these 
equations and compared with experimental values in figures 1 6 
and 17. The .general trend of the experimental points for’ t « 3 d 
is -higher than that for t »-'12 0 . This is in agreement with 
theoretical prediction. Either of these curves may be used for ‘ . 

design purposes. . 

For research purposes equations (79) and (82) may "be written 
in the following forms! . 



In each case the first definition of contains quantities which 
may "be measured experimentally. The second definition is a function 
of r Q only. The i]/ functions are universal functions of r 0 
being in no way affected by variations in t or p . These \{/ 
functions have been plotted and- compared with experiment in' 
figures 18 and 19 ^ These .functions are not impendent upon the 
properties of the float or the definition of the additional water 
mass. 
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' CONCIDSIONS 


The landing impact of a seaplane may he analyzed hy means of 
Newton's lav of motion and elementary hydrodynamic theory, providing 
a number of restrictive assumptions are made. Some of these 
assumptions ai'e as follows i 

1. The water surface is smooth* 

2. The seaplane is a rigid body. 

■ 3 . The trim angle remains constant during the impact. 

1*. The wedge bottom of the hull is not flared. 

5* The finite width of the hull has no effect on the motion. 

6 . The wing lift equals the weight of the seaplane. 

7. Buoyancy and viscosity forces may be neglected. 

The three-dimensional apparent mass of the water may be taken as 
a direct analog of the two-dimensional concept of elementary hydro- 
dynamics. On this be$is a solution for maximum acceleration may 
be obtained which snows good agreement with experimental results. 
Agreement between theory and experiment in regard to displacements 
is more difficult to obtain and calls for further improvements in 
the analysis through theoretical and experimental research. The 
TneTrimum acceleration 1 is found to be proportional to the square of 
the initial velocity and inversely proportional to the draft at 
the instant of maximum acceleration. The constant of proportionality 
is given by the theoretical analysis and is found to be independent 
of the float properties or the definition of the additional water 
mass. 
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APPENDIX 


In the preparation of the graphical Illustrations of .the paper 
a number of tables of functions were computed. Some of the tables 
are included in this appendix in order to provide assistance in 
further . researah studies and for design use of the graphs of 
accelerations. The following function -is given in table' I. 

(f = log (1 + r) + -. A— (100) 

1 + r 

In table IT are given the values of r m and as obtained by 

simultaneous solution of equations (86) and (87) » The values of r m 
were read from the intersections of graphs and cannot be considered 
as accurate in the last significant figure. The values of were 

not read from the graphical intersections but were computed from . 
equation (87) using the values of r ffi . 

In table III the values of u n are given as obtained by 
solving equation (95) • In table IT values of are given as 

defined by equation (89) . Table T gives values of the functions 
\J /1 and ^2 as defined by equations (98) and (99).' 
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TABLE I 


r 

4 

r 

4 

0 

1.00000 

2.1 

1.45398 

' .1 

l.oo44o 

2.2 

1.47565 

.2 

1.01565 

2.3 

1.49695 

.3 

1.03159 

2.4 

1.51790 

.4 

1.05076 

2.5 

1.53847 

.5 

1.07214 

2.6 

1.55871 

.6 

1.09500 

2.7 

1.57860 

.7 

1.11887 

2.8 

1.59816 

.8 

1.14334 

2.9 

1.61739 

• 9 

1.16817 

3.0 

1.63629 

1.0 

1.19315 

3.1 

I.65489 

lrl 

1.21813 

3.2 

1.67318 

1.2 

1.24301 

3.3 

1.69118 

1.3 

1.26769 

3.4 

1.70887 

1,4 

1.29214 

3.5 

I.72630 

1.5 

1.31629 

3-6 

1.74345 

1.6 

1.34013 

3.7 

1.76033 

1.7 

1.36362 

3.8 

1.77695 

1.8 ' 

1,38676 

3.9 

1.79332 

1.9 

1.40954 

4.0 

1.80944 

2.0 

1.43194 
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TABIE II 


r 

0 

r m 

m 

0 

0 

0 ■ 

.2 

0 

tr\ 

0 

* 

^0157 

.4. 

• .149 

.0423 

.6 

. . .277 

• .0698 

.8 

• ' .421 

.0941 

1.0 

•361 ■ 

.1130 

1.2 

-m 

.1294 

1.4 

.660 

.1431 

-1.6 

1*€12 

.1547 

1.8. ■■ 

1.166 

.1647 

2.0 - 

. 1.318. 

} 

.1731 

2.2 

1.471 

i 

.1805 

2.4 

1.626 

.1871 

2.6 

1.782 

. .1929 

2.8 

1.937 

.1981 

3-0 

2.090 

.2026 

3.2 

. 2.244 

.2067 

3*4 

i 

2.400 

.2105 

3.6 

2.554 

.2139 

3-8 

2.710 • 

.2171 

4.6 

2.865 

.2199 
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TABLE 17 

.... VALUES QF" oc£ ' 
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TABLE Y 
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Figure I.— Seaplane aV )ns-Fan4- of- En+ry 
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Ficjure DispiacernervVs and R-eac+ion 

Durin9 Innpact Period 
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(c) NATIONAL ADVISORY 
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Fiqure 3 — Reac+ion Distribution 
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Z. NATIONAL ADVISORY 

COMMITTEE FOB AERONAUTICS 


Figure 4-.— ELlemervfaru Concept 
of Additional Mass 
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Figure 5 — Variations with Displ£Lcemen+ 
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Figure II— Displacement at Ins+ant of Maximum Acceleration 
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Fiqure 14". ~ U.on^'i'Vuciina^ Dis'tribU'Vion 

of Reac+ion 
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Fiqur'e 15. — Maximum Drakf+ 
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Figure 16.— Rak+io erf x-rJ+o 
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